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EFFECT  OF  VARYING  AIR  DENSITY  ON 
THE  NONLINEAR  PITCHING  AND  YAWING  MOTION  OF  A  SYMMETRIC  MISSILE 

ABSTRACT 

Although  the  effect  of  an  exponentially  varying  air  density  on 
the  angxilar  motion  of  a  re-entering  missile  has  been  intensively 
studied  for  linear  aerodynamic  moments,  this  effect  for  nonlinear 
moments  has  only  been  investigated  by  a  quasi-linear  technique  and 
this  work  was  limited  to  planar  motion.  In  this  report,  various 
modifications  of  the  quasi-linear  analysis  for  a  cubic  static  moment 
are  compared  with  the  exact  solutions  and  the  best  of  these  used  to 
determine  the  general  effect  of  varying  density  on  combined  pitching 
and  yawing  motion.  Next  the  perturbation  method  which  makes  use  of 
the  exact  solution  for  a  cubic  static  moment  is  considered  and  the 
necessary  modifications  introduced  by  the  varying  density  are  made. 

The  predicted  density- gradient- induced  damping  for  planar  motion  is 
determined  for  both  theories  and  compared  with  the  numerical  integration 
of  the  exact  differential  equation.  The  perturbation  prediction  is 
found  to  be  superior. 
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I  _i  >»  ,  p 

twice  the  total  energy,  |  +V(&) 


C 


2 


D 

E(k) 

H 


K(k) 

K,,  Kg 


twice  the  trajectory  component  of  the  angulsir 
momentum  of  the  pitching  and  yawing  motion, 

Cg  =  i(I*l  -  i  I) 

drag  coefficient 
lift  coefficient 

static  moment  coefficient 

damping  moment  coefficients 


Magnus  moment  coefficient 


coefficients  of  the  transverse  components  of  the 
aerodynamic  moment 

*  2 
non-conservative  part  of  C,,  ,  C  =  C,,  -  c  -  C-B 

M  '  M  o  2 

a  a 

,  2»2ft-2 

to  oTM 
o 

complete  elliptic  Integral  of  the  second  kind 

•psi  r  r,  r,  -2 


2m 


+  )' 
% 


axial  moment  of  inertia 


transverse  moments  of  inertia 


complete  elliptic  integral  of  the  first  kind 
amplitudes  of  two  modes  of  linear  oscillation 
modulus  of  the  elliptic  Integral 
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u,  V,  w 
V 

X,  y 
z 

a 

P 

7 

6 

^2 

e 

e 

I 

r 

t 

p 

PQ 


components  of  velocity 

AfS  2~  i" 

magnitude  of  velocity,  V  =  \^u  +  v  +  w 

A 

the  potential  function  associated  \/ith  M 

squared  amplitudes  of  the  two  modes  of  linear  oscillation 

altitude  (ft.) 

angle  of  attack 

angle  of  sideslip 

cosine  of  total  angle  of  attack 

IM 

minimum  value  of  6 
maximum  value  of  8 

angle  the  flight  path  makes  with  respect  to  the  vertical 
argument  f 

aerodynamic  damping  doefficients 

V  ^  iw 

V 


air  density 

sea- level  air  density 

coefficient  of  exponential  density  variation 
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a 

aicosO 

phase  angle  of  the  jth  mode 

K-h 

h 

a? 

Subscript 

a 

average  value  over  cycle  of  motion 

Superscripts 

t 

derivative  with  respect  to  arclength 

- 

complex  conjugate 

quantity  related  to  non-rotating  coordinate  system 


quantity  related  to  coordinate  system  which  is  rotating 
p 

with  velocitv  •= 
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1.  DlTRODUCTIOx'T 


The  influence  of  an  exponentially  varying  air  density  on  the 

planar  pitching  of  a  re-entering  missile  with  linear  aerodynamic 

1  2 

moments  has  heen  studied  "by  a  number  of  authors.  *  This  work  was 

extended  to  combined  pitching  and  yawing  motion  of  a  spinning  missile 

5  h  , . 

by  Leon  for  no  aerodynamic  damping  and  "by  GLU'ber  for  linear 

5 

aerodynamic  damping.  Recently,  Coakley,  Lail.one,  ana  :'iiss  have  made 

(• 

use  of  a  quasi-linear  technique  to  describe  the  planar  motion  of  a 
missile  with  a  cubic  static  moment  flying  throuijh  an  exponential 
atmosphere. 

In  this  report  we  will  study  the  influence  of  such  density 
variations  on  the  general  combined  pitching  and  yawing  motion  of  a 
missile  acted  on  by  nonlinear  moments.  The  nonlinear  analysis  will 
make  use  of  a  perturbation  technique^  which  is  more  accurate  than  the 
quasi-linear  analysis  employed  by  Coakley,  Laitone,  and  Mass.  The 
quasi-llneeir  analysis  will,  however,  be  derived  for  combined  pitching 
and  yawing  motion  for  comparison  with  the  more  exact  method.  Although 
the  primary  objective  is  the  influence  of  varying  density,  the 
development  will  be  so  formulated  that  the  effect  of  Mach  number  or 
Reynold* s  number  variation  of  the  aerodynamic  coefficients  themselves 
may  be  studied. 


2.  QUASI-UNEAR  SOLUTION  FOR  CONSTANP  DENSITY  AND  NO  DAMPING 

In  this  section  we  will  consider  various  raodif icatlons  of  the 
quasi-linear  technique  for  no  aerodynamic  damping  or  drag,  constant 
density,  and  cubic  static  moment.  For  this  case  the  aerodynamic 
moment  may  be  written  in  the  form: 

My  +  iM^  =  -  (i/2)pV^S/(c^  +  (1) 

where  6  =  =  p+ia 

8  =  I  I  I  and  the  other  symbols 
are  defined  in  the  List  of  Symbols. 


In  Equation  (l)  the  complex  moment  and  complex  angle  are  in  a  missile- 
fixed  coordinate  system.  .If  we  derive  the  differential  equation  for 
the  pitching  and  yawing  motion  for  tills  moment*^  in  a  non-rotating 
coordinate  system,  it  would  have  the  form^ 


cosine  of  total  angle  of  attack 


6"  -  (  7  +  ^7  5*  -  (M^  +  =  0 

where  6  =  6  exp  ^  dsj 
7  =  (1  -  8^)^/^ 

P  =  5l  Pi 

-•■k)  •• 

¥)  [ 


Cg  -  (1/2)C^ 


] 


and 


(2) 


dt  is  the  Independent  variable. 


*  Since  the  static  moment  coefficient  is  multiplied  by  7,  a  cubic  moment 
produces  higher  powers  in  the  equation  of  motion.  Only  the  cubic 
powers  are  retained  in  Equation  (2). 
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For  moderate  geometrical  angles  y*  =  0.  For  constant  P,  a  transformation 
to  a  coordinate  system  which  has  an  angular  rate  of  P/2  reduces,  this 
equation  to  a  simple  form. 


I"  -  5  (1  +  m_8^)  t  =  0 


(5) 


where 


I  i(l/2)Ps 


A 

M  =  M 
o  o 


-  i^A 


“2  =  V^o  • 

For  a  linear  moment  and  negative  |  moves  along  an  ellipse 
with  equation 


A 

6  =  Kj^e  +  K^e 


(M 


where  K,  are  constants 
J 


.  +  S'.s 
Jo  ^J 


and 


?]■-  fA  . 

The  semi-major  axis  is  and  the  semi-minor  axis  is  j  -  K^j 

When  the  spin  is  zerp  the  negative  requirement  is  that  for  static 
stability.  In  the  case  of  a  spinning  missile,  this  requirement  is 
essentially  that  for  gyroscopic  stability 


<1 


(5) 


where  s  = 

g 


p2 

Tm 
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The  motion  in  the  non-rotatins  ii'arne  of  reference  becomea  epicyclic. 

The  frequencies  in  these  c  .jvdiuul.e;'. ;  5^.^,  fev  a,  sally  stable  missile 

(M^  <  O)  ^6  opposite  in  sItu  aiui  for  a  itacie.alij  u.i  d.tbie  missile  they 
have  the  same  sign  as  that  of  f. 

The  quasi-linear  method  assumes  a  solution  f  -a  iv;  ■  iju  (  )  ^  ■  the 
form  of  Equation  (4)  with  frequencies  which  depend  ;  i,.; .  :  r 

actual  calculation  of  this  dependence  malcea  use  of  algi  bi*a  used  in  the 
method  of  variation  of  parameters.  If  Equation  (>)  ..rlLben  in  the  form 

t'  -  (ij  (6) 

Equation  (4)  is  the  solution  for  =  0.  We  assume  the  solution  of 
Equation  (6)  to  be 


4  A  '  4 

?  =  Kj^e  +  K^e 


(7) 


Then 


■'j  ■ 


coo  ^ 

where  % 

V  f— r-  r  ‘51,  41  .  .  4 

1  =  i  \1  -  iK^e  ^  -  KgO  +  (K^  +  ^ 


(0) 


/  •  * 

+  (143+ 


(9) 


I  I  I  I 

Let  (Kj^  +  itj^K^)e  +  (K^  v  "•  =  0  (lO) 


l4 


and  differentiate  again 

6  =%\K^e  J 


+  i  \j  -  |^(K^  +  i>lf3_Kj^)e  -  (K^  +  itgKg)® 

Equations  (7  -  U)  are  now  substituted  in  Equation  (6). 


Equation  (12)  and  a  similar  equation  for  the  other  mode  are  exact  but 
rendered  quite  complicated  by  the  presence  of  The  basic  assumption 
that  damping  and  frequency  shift  over  a  cycle  of  ^  is  small  has  to  be 
made.  If  this  is  the  case  Equation  (12)  may  be  averaged  over  a  cycle 
of  ^  with  the  resvilt 

k'  .  %(!^  t  2Il|)  \1-  M 

iq  *  ‘h  =  ‘  — 2-  ■  '*  °  («1 

Thus  the  quasi  linear  solution  has  no’  damping  but  does  have  frequencies 


%  ■  -  \l-'l  [i  *  ^  (2!^  »  4)j 


This  method  can  easily  be  extended  to  treat  both  linear  and  nonlinear 

7 

damping  . 

Equations  (l4-15)  can  be  substantially  improved  in  accuracy  if 
the  average  value  of  the  coefficient  of  5  is  used  on  the  left  side  of 
Equation  (6)  instead  of  its  linear  value.  Since  the  average  value 
of  6^  is  +  K^,  Equation  (6)  would  become 

-  M,  ji  +  t=  ^  -  0 
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From  this  equation,  quasi-ltnear  vaiiies  of  ^ .  may  he  derived  in  a 

J 

manner  similar  to  the  above . 


This  improved  quasi  linear  technique  was  used  with  some  success  in 
Reference  6.  One  clear  advantage  of  Equations  (17-I8)  is  that  they 
may  he  applied  to  periodic  motion  of  a  missile  which  is  gyroscopically 
\uistahle  for  small  amplitudes  0,  nig^O).  In  this  case  infinite 

frequencies  are  predicted  for  ^  ^  ~  As  we  will  see  the  correct 

answers  are  hounded.  To  avoid  this  difficulty  and  obtain  sli^tly  more 
accurate  estimates  for  other  cases,  the  substitution  method  of  Reference  8 
can  he  used.  This  method  will  he  developed  as  another  variant  of  quasi- 
linear  method. 


We  assume  a  solution  of  the  form 

.  % 

where  K  and  ^  are  fvuictions  of  s. 
J  J 


Then, 


Let 


Tills  equation  is  essentially  the  neglection  of  damping  in  comparison 
with  frequency.  We  now  differentiate  again. 


1"  =  1  ^  ^ 

Ati  Ati  l-^gl 

+  w  +4"^^  J 

Equations  (19  -  22)  can  now  be  substituted  in  Equations  (3)  and  the 


result  manipulated  into  the  form 


Q'f  -  !(?;  -  ?2)  I  -  -  8,,  -  ^  e-l-?)  ^ 

-  \i%f  *  K  ■  €]  J  *  € 

For  linear  moments  with  constant  coefficients  the  right  side  of 
Equation  (23)  is  constant.  When  the  moment  is  nonlinear  or  the 
coefficients  vary,  approximations  for  frequency  and  damping  may  be 
obtained  by  averaging  the  right  side  over  a  cycle  of 

K  =  ^  [  1  +  "2(1^  +  24’]  < 


Similar  relations  apply  for  the  other  mode 


When  or  trig  are  functions  of  s.  Equations  (25)  and  (27)  allow  us  to 

compute  the  effect  of  this  on  the  damping.  This  will  be  done  in 

/> 

Section  3-  For  constant  and  damping  is  zero  and  the  frequencies 
are  given  by  Equations  (2k)  and  (26). 
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With  these  three  sets  of  values  for  the  f rciiueiicies  derived  it  is 
quite  important  to  determine  their  i*elative  accuracy.  To  do  this  we 
will  make  use  of  the  exact  solution  for  8.  Four-  fundamentally  different 
variations  of  static  moment  may  be  expressed  by  Equation  (p)  wiien 
different  algebraic  signs  are  assigned  to  the  coefficient  and 
Onl^  three  of  these  moments,  however,  can  cause  periodic  motion.  These 
three  are  illustrated  in  Figure  1  and  may  be  identified  in  the  following 
way: 

(a)  Stable  at  small  angles;  more  stable  at  larger 
angles  ">2  >  0) 

(^)  Stable  at  small  angles;  less  stable  at  larger 
angles  <0,  nig  <;0) 

(c)  Unstable  at  small  angles;  stable  at  larger  angles 

(^o  >  0^  “2  <0) 

The  periodic  solutions  of  Equation  (5)  are  derived  in  References  6 
and  9  and  are  sumarized  below 

type  (a)  moment 

5^  =  82  -  (82  -  8^)  sn^(ffis,  k)  (28) 

type  (b)  moment 

5^  =  8^  +  (8|  -  6^)  k)  (29) 

4 

2 

CO  >  0 

type  (c)  moment 

8^  =  62  -  (82  -  8^)  sn^((os,  k)  (30) 
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where  is  minimum  value  of  B 


Bg  is  maximum  value  of  8 


The  ineq.ualities  associated  with  Equations  (29-30)  are  quite 
important  in  themselves.  Three  interesting  observations  may  be  made: 

1.  For  type  (a)  static  momert,  periodic  motion  of 
any  amplitude  is  possible. 

2.  For  type  (b)  static  moment,  symmetric  planar 

periodic  motions  (P  =  0,  =  0)  are  possible  for  all 

amplitudes  for  which  the  moment  does  not  change  sign; 
circular  motion  (8^  =  Bg),  however,  is  possible  only  when 
the  nonlinear  part  of  the  moment  is  not  greater  than  two- 
thirds  of  the  linear  part. 

3.  For  type  (c)  static  moment,  possible  periodic 

2 

motions  are  those  for  which  the  median  value  of  8 
yields  a  stable  moment. 


The  veiriables  8^  and  8g  may  be  approximately  related  to  the  modal 
amplitudes  and  Kg  by  the  equations 

®1  =  ^^1  ■  (51) 

=  (^1  (52) 
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The  pci'iod  of  the  elliptic  cine  function  in  L’quid, ion:;  (28  -  JO) 
corresponds  to  half  tlie  perici  of  '■(>  i:i  Eiiuatio/i  (o)  uiucc  thCib  equation 
can  be  put  in  the  form 

.2  /„2 


^  =  ^2  ■  ^^2  “  sin"-(i^/2) 

=  6^  +  (Sg  -  sj)  sin^  (tJ  +  «)/2 


f 


m 

where  K(k)  is  the  complete  elliptic  intecral  of  the  first 
kind.  (The  period  of  sm  (ws,k)  is  hK/u. ) 

The  vaxious  approximate  values  of  %  may  now  be  compare!  with  the 
exact  value  of  Equation  (34).  In  Figures  2-3  this  is  done  for  the  three 
different  moment  types  and  both  planar  and  circular  motion.  As  can  be 
seen  from  these  plots  of  {  versus  m  =  tne  unaltered 

quasi-linear  estimate  (QL)  is  only  good  for  small  amplitudes  while  both 
the  improved  quasi-linear  (IQL)  and  the  substitution  methods  are  quite 
good  for  the  reinge  of  m  considered  with  the  exception  of  the  interval 
-  1  ^  m  1/2.  Since  is  positive  for  a  type  (c)  moment,  a  quasi- 
linear  value  of  ^  can  not  be  computed.  Although  substitution  method  is 
only  slightly  better  than  the  improved  quasi-linear,  its  potentiality  of 
describing  the  effect  of  varying  coefficients  and  nig  make  it  more 
valuable  than  the  Improved  quasi-linear.  The  good  agreement  of  both 
methods  encourages  us  to  introduce  varying  coefficients  as  veil  as 
nonlinear  damping. 


5.  QUASI-miEAR  SOU/riON  FOR  VARYING  DENSITY 


The  coefficients  in  the  differential  equation  for  a  jnissile*s 
pitching  and  yawing  motion  may  he  functions  of  the  independent  variable 
for  a  nujid)er  of  reasons:  varying  air  density,  varying  missile  mass 
and/or  moments  of  inertia,  varying  aerodynamic  coefficients  as  a  result 
of  their  dependence  on  Mach  ntimber  or  Reynold* s  number.  At  the  present 
the  greatest  interest  lies  in  the  influence  of  varying  density  on  a 
missile  leaving  or  entering  the  earth’s  atmosphere.  Althou^  most  of 
the  equations  of  this  report  will  apply  to  any  cause  of  varying 
coefficients,  the  effect  of  density  will  be  our  primary  objective. 


The  actual  density  variation  can  be  reasonably  well  approximated* 
by  an  exponential  up  to  altitudes  of  300, OCX)  ft. 


(35) 


=  sea  level  density 
z  =  altitude  in  feet  and 
®  “  22,000  ft 

If  9(3)  is  the  angle  the  flight  path  makes  with  respect  to  the  verticaQ., 
z  can  be  related  to  our  independent  variable  s  by  the  equation: 

s 

z*  =  -  ^  i  cos  0(sj^)  dSj^  (36) 

o 

where  0  <  9  <  90^  entering  the  atmosphere 
90°<  0  <  180*^  leaving  the  atmosphere 


Domraett  '  has  approximated  the  ARDC  model  atmosphere  by  a  set  of  four 
exponentials.  This  more  accurate  description  could  be  used  in  the 
theory  of  this  report  and  is  described  in  the  appendix. 


Eq.uatijns  (55“36)  can  now  be  used  to  obtain  derivatives  of  and 


M*  or  M 
o  _  _ o 

ft  ft„ 

o  o 


nig  ^ 


"fe  A 


a 

1  -  s 


o  ^ 
—  =  -  o 


s 

g 

1  -  s 

g 


where  a" =  at  cos  o  =  —  • 

P 

Note  the  simple  form  these  derivatives  assime  for  zero  spin,  (s  =  0) 

s 

With  the  introduction  of  varying  ft^  and  we  can  now  return  to 
Equations  (24-27)  which  were  derived  by  the  quasi-linear  theory  and 
consider  the  effect  of  this  variation.  In  order  to  increase  the 
generality  of  the  results  both  linear  and  nonlinear  aerodynamic  damping 
will  be  Introduced.  Thus  Equation  (3)  will  be  replaced  by 

i"  -  A,j(l  +  BgS®)  1  •  -  Bt'  +[m*  +  1P(T  -  |)lt  (59) 

vhere  H  -  ('m  *  ^  '=«.>! 

a  1  a  J 


_  psi 


m  pSi 
T  =  ^  7 


-2  * 


p  2m  t  M 


J  :z 


is  transverse  radius  of  gyration 


^a  axial,  radius  of  gyration  and 

=*-  0„_^.o„.=/.cV,  (eV) 


It  shoiild  be  emphasized  that  the  aerodynamic  coefficients  defined 

in  Eq.uation  (39)  coefficients  and  not  derivatives.  They  may  be 

functions  of  6^  and  (6^)  .  C*  is  the  nonpotential  part  of  Cj-  .  Since 

a 

this  is  a  rather  strange  quantity,  we  will  consider  the  value  of  it  and 
the  other  quantities  in  Equation  (39)  for  a  nonspirming  missile  with 
cubic  moments  8uid  forces.  (Example  3  of  Refei*ence  6).  The  force  and 
moment  were  assumed  to  lave  the  form 


°D  =  ®o 


ffl  n 


•a*" 

0.  +  10„  =  -  1  +  Cg**  +  * 


<=M  *  V  '  ^  -  'U>^ 

q  a 

C*  =  (o,,)  (6^)' 


u*  - 

“  =2ir 


[(a  -  e  -  k”^d  )  +  (a„  - 

o  o  t  o^  ^  2 


If  the  small  effect  of  aerodynamic  damping  on  frequency  is 
neglected,  the  substitution  quasi-linear  method  yields  the  following 
relations . 


'4-  -\1  -  t  *  "2K  + 


Xl 

t 

h 


ji. 


K-% 


%-K 


+  X, 


m 

m 

(50) 

(51) 


where 


2«{^i  -  Sa) 


^2  ^ 


] 


cos?  7 


-  M*  sin  ^  -  P(T  -  |)  ^ 

Equations  (48-49)  may  be  differentiated  and  solved  simultaneously 
for'$(^/(0^  -  {^g)  and  ^g/ (02  -  0ji). 

K  ^  ^ 

=  ^  aqi  +  2Xj^a^  +  2Xga^^  +  —  (a^  +  (52) 

M  ^ 

o 


-  S!, 
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Al 


®21 


+  ^  (a. 


‘22 


+  agj)  (53) 


where  a.,  are  defined  in  Table  I. 

jk 

Equations  (52-53)  can  be  used  with  Equations  (50-51)  to  predict 
the  variations  in  aiiqjlitude  but  the  eilgebraic  expressions  are  quite 
cumbersome.  For  two  special  cases  of  some  importance  considerable 
simplification  is  possible.  The  first  case  is  that  of  a  linear*  moment* 
(iQg  =  O)  and  the  second  is  motion  through  zero  amplitude 
(®1  **■  11“®“^  moment  Equations  (50-51)  reduce  to 

the  very  siiiq)le  form 


*1 


At 

M 

o 


a 


-  (I/I*)  ^  =  »-i  -  Mi  -y 


(54) 


i  -  Xj  -  (V4)l^)-  4,  -  “ 

.  Qi 


(55) 


According  to  Equations  (5^-55)  the  aerodynamic  danqxing  (\j)  of 
a  reentering  missile  will  grow  from  zero  to  sea  level  values,  while  the 
Influence  of  the  density  gradient  (a)  depends  on  the  stability  factor. 

The  motion  of  a  statically  stable  nonspinning  missile  will  be  damped 

« 

by  this  term.  If  it  possesses  constant  nonzero  spin*  this  density 
gradient  damping  will  grow  from  zero  as  s  varies  from  minus  infinity 
to  its  sea  level  valuel  A  staticsilly  unstable  missile  with  constant 
spin  win  have  density-gradient-induced  undein^xing  which  grows  from  zero 
as  Sg  decreases  from  plus  infinity. 


McShane,  Kelley,  and  Reno^  consider  the  effect  of  varying  spin 
by  means  of  the  WKB  method. 
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This  strange  state  of  affairs  may  be  quickly  resolved.  The 


quasl-linear  technique  assumes  that  changes  over  a  cycle  are  reasonably 
small,  l.e.,  is  small.  This  condition  determines  a 


is  small.  This  condition  determines  a 


maximum  altitude  for  which  the  theory  is  applicable.  For  some  missiles 
a  maximum  altitude  of  300,000  ft.  is  appropriate.  At  Lhis  altitude  the 
effect  of  increasing  spin  can  be  determined  and  as  caia  be  seen  from 
Equations  (5*^-55)  the  damping  for  a  statically  stable  missile  and  the 
undamplng  for  a  statically  unstable  missile  are  both  reduced. 


When  the  pitching  and  yawing- motion  goes  throu^  zero  amplitude 

A,»  /b»  ^2 

=  -  ^2  ^  ^  ~  relations  are  very  helpful  in 

reducing  Equations  (52-53)  to  a  reasonable  degree  of  complexity. 
Equations  (50-51),  then,  become 


If  ^  will  not  remain  equal  to  K^.  Equations  (56-57), 

therefore,  are  valid  only  when  is  nearly  equal  to  or  the  Tnin-ImiiTTi 
6  is  near  zero. 
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For  plansir  motion,  spin  is  zero  and  the  aerodynamic  damping 
rates  are  equal.  Equations  (56-57)  col3.apae  into  a  single  equation. 


(^-1) 

Coakley,  Laitone  and  Mass^  derived  Equation  (58)  for  1  =  0  in  a  somewhat 
different  manner.  The  pole  at  m  =  -  8/9  has  rather  strange  consequences. 
According  to  Eqviation  (58)  the  angular  motion  of  a  missile  with  a  type  (b) 
moment,  ascending  in  the  earth's  atmosphere,  will  grow  in  amplitude  until 
M 

5-  =  -  (8/9)  •  This  upper  bound  does  not  exist  for  'J'  =  0  since 

d.  Mn 

Equation  (29)  has  solutions  for  -  1  <  m  <  0.  This  difficulty  will  be 
resolved  in  the  next  section. 
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4.  PERTURMTION  SOLUTION  FOR  VARYING  DENSITY 

In  Reference  6,  the  nonlinear  damping  is  treated  by  perturbing 
the  exact  solution  of  cubic  static  moment  equation.  As  is  shown  in 
Section  2  this  solution  Involves  an  elliptic  function.  The  perturbation 
method  makes  use  of  two  quasi-constants  of  the  motion,  the  energy  and 
angular  momentum. 

=1'  r  +v(8^)  =  (5')^  +  (80*)^  +v(8^)  (59) 


c 


2 


Kf 


A*  K.  2a* 

I  6)  =  28‘'9 


(60) 


where  V(8  )  =  -  11^(8 


Eqviation  (39)  ®ay  now  be  rewritten  in  terms  of  derivatives  of  these 
quasi-constants 


=  *^2H(Cj^  -  ^)  +  (M*)(8^)'  +  P(T  -  |)  Cg 


C*  =  -  HCg  +  2P(t  -  |)  8^ 


(62) 


A  *  2  /  2v  * 

H,  V,  M  and  T  are  functions  of  8  and  (8  )  .  These  quantities  are 

« 

periodic  functions  given  by  Equations  (28-30).  The  right  sides  of 

Equations  (61-62)  may  be  averaged  over  the  period,  P  ,  to  yield  functions 
2  2 

of  Cg,  8^,  and  82.  Generalized  modal  amplitudes  may  now  be 
introduced. 


8^  =  (yr-x  -  yrjf  (63) 

=  (  vir+  vyf  (64) 
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Tha  quasi-constants  may  be  expressed  in  terms  of  the  generalized  modal 
amplitudes 

CjL  =  -  [[2(x  +  y)  +  ^  (3x‘^  +  lOxy  +  5y^)J  (65) 


Cg  =  2(x  -  y)  ^  -  Mq(1  +  iDgCx  +  y))  (66) 

Equations  (65-66)  can  be  differentiated  and  solved  for  x*  and  y*. 

[2  +  nu(x  +  3y)]  (c*)*  +  (2  +  nu(5x  +  3yj]  N1-  (1  +  nu(x  +  y))  (C*)^ 

-  2D1? 


[2  +  n^(3x  +  y||  (C^)*  -  (2  +  b^(3x  +  5y3\|-  ^^2^* 

’’’  -  2I)fl 


where  _  .A-S  2-12 
D  =  4M_  m  <a 


=  j2  +  d^(3x  +  2  +  3y2l  [2  +  n^(3x  -  2  +  3y)3 

=1)*  -  (^) 


(o’)* .  o'  -i  r°  ^  c  -  I  i) 

2  2  2  "aj  2  \[-  (J^(l  +  m^(x  +  y)) 

t  ^  t  # 

(Cj^)  and  (Cg)  may  be  evaluated  by  the  use  of  Equations  (61-62) 


'  i  po  ^  "2]  c 


30 


(< 

+  P(T  -  f)  ^  +  SJ^  ||\  p  -  2{x  +  y)j- 
«4)*  .  ^  j^H  0^  *  2P(T  -  |)  eO  ^ 


l(x  -  y) 


(1  +  nv,(x  + 


idiere 


n.-i-  W] 


For  a  Donspinnlng  missile  in  on  exponential  atmosphere  is  zero  and 
the  effect  of  the  density  gradient  is  to  replace  H  by  H  +  ^.  In  fact 
an  inspection  of  the  derivation  of  proceeding  equations  indicates  that 
this  replacement  of  H  by  H  +  ^  for  a  nonspinning  missile  in  an 
exponential  atmosphere  is  true  for  any  ±.e.,  any  nonlinear  static 
moment.  Here  again  it  must  be  eBq)hasized  that  the  results  of  the  theory 
are  valid  only  when  the  density- gradient- induced  damping  is  small  over  a 
cycle  of  the  basic  periodic  motion.  The  averages  of  for  circular 
motion  are  given  in  Reference  6  and  for  planar  motion  in  Reference  11. 
For  planar  motion  (x  =  y  =  (1/4)  6^,  =  O)  and  a  constant  H, 

Equations  (67-68)  collapse  to 


(2  +  ni  -  2/^2  -  mAj^) 


(1  +  m) 


[s-r]  • 


For  a  type  (b)  moment 


*2'  V 


b  -  "a] 

[2(1 + 1 


A4  =  (V5)  kf  [2(1  +  kjfjAg  -  ij 


-  K(kp) 


conplete  elliptic  integral  of  the  first  kind 


E  =  E(k  ) 
P  P 


comp3.ete  elliptic  integral  of  the  second  kind 


k^=  -  ° 

p  2  +  m 


modulus  for  planar  motion 


The  coefficient  of  •jj-  in  Equation  (7I)  has  a  pole  at  m  =  -  1. 

Since  this  is  the  upper  bound  for  Sg  io  Equation  (29),  this  is  much 

more  reasonable  than  the  pole  at  m  =  -  8/9  in  Equation  (5S)*  The 
‘cT 

coefficient  of  in  these  two  Equations  is  plotted  versus  m  in 
Figure  4.  Equation  (3)  was  numerically  integrated  for  planar  motion 


and  an  exponential  density  and  the  logarithmic  derivaWve  of  6g 
calculated.  Numerical  values  of  the  coefficient  of  ^  could  thus  be 


determined  and  are  plotted  in  Figure  4.  The  agreement  with  Equation  (7I) 
is  very  gratifying. 


5.  SUMMARY 


1.  The  improved  quasi-linear  and  substitution  methods  for  a 
cubic  static  moment  predict  frequencies  to  very  good  accuracy. 

2.  These  methods  do  not  predict  the  effect  of  varying  density 
to  quite  the  same  accuracy. 

3.  The  effect  of  density  gradient  on  the  damping  of  missiles 
with  linear  static  moment  is  to  add  3/^1  -  s^)  to  each  dating 
exponent. 

4.  The  effect  of  density  gradient  on  the  damping  of  a  nonspinning 
missile  with  a  nonlinear  static  moment  is  to  add  3/2  to  H.  This 
modification  allows  the  use  of  the  results  of  the  perturbation  theory. 

CHARLES  H.  MURPHY  /  ^ 
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APHMDIX 

A  QUADRI-EXrONEfflllAL  A'tMOSPIffiKE 

In  Reference  12,  the  ARDC  model  atmosphere  is  approximated  by  a 
set  of  four  exponentials  of  the  form 


“g 

where  =  I.0l<-1  x  lO”^  oK/ft 

X\ 


z  =  0  and  the  other  constants 
o 

are  given  in  the  Table. 

The  Vf  lues  of  a  given  in  the  Table  could  be  used  in  the  theory  of  this 
report  at  the  appropriate  altitudes  instead  of  1/22,000  ft. 


TABT.K’^ 


"i 

.V“i 

Pi 

■2 

'^i 

ft. 

ft. 

^  slu2s/ft. 

°K 

35,000 

30,800 

2.377  X  10"^ 

288.16 

11jO,000 

21,000 

4.034  X  10"^ 

218.91 

2lt0,000 

26,900 

9.471  X  10"^ 

279.68 

300,000 

18,600 

5.099  X  10"^ 

193.41 

*  These  values  are  rounded  .from  those  of  Reference  12. 
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